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Abstract 

We find approximate analytical presentation of the solutions 9(ri,r2,rx2) °f Schrodinger equa- 
tion for two-electron system bound by the nucleus, in the space region 7*1,2 = and r\% = that 
are of great importance for a number of physical processes. The presentation is based on the 

in ' 

well known behavior of VE^ri, r^i t\%) near the singular triple coalescence point. The approximate 

o : 
• i— i 

ion of hydrogen H . The functions are shown to determine a natural basis for the expansion of 



functions are compared to the locally correct ones obtained earlier by the Correlation Function Hy- 
perspherical Harmonic (CFHH) method for helium atom, light heliumlike ions and for the negative 



X 
S3 



CFHH functions in the considered space region. We demonstrate, how these approximate functions 
simplify the calculations of the high energy ionization processes. 
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I. INTRODUCTION 



The ground states of the two-electron systems bound by the nucleus are described by 
the wave functions, depending on three variables. These can be the distances between the 
electrons and the nucleus r 12 and the interelectron distance r 12 . Here we shall find analytical 
expressions, which would approximate the solutions of the Schrodinger equation \?(ri, r 2 , ri 2 ) 
in the special cases r 12 = and r 12 = 0: 

F(R) = (0, R, R); <$>(R) = ^(i?, R, 0) , (1) 

We shall consider the ground states of helium atom and of the light heliumlike ions, including 
also the negative ion of hydrogen H~ . In this paper we shall treat the ground states only 
Thus, the total spin of the two-electron system is equal to zero. 

Note that this problem differs strongly from the traditional problem of approximating the 
total wave function ^f(ri, r 2 , r 12 ) There are numerous wave functions of this kind, with 
the approximate functions being usually certain combinations of exponentials and polyno- 
mials, while a set of fitting parameters is found by minimization of the energy functional. 
Thus, the quality of such functions is determined by the accuracy of reprodusing the binding 
energy value. Since the averaged value of the Hamiltonian is determined by the distances, 
which are of the order of the size of the atom, such functions provide very good approxima- 
tion at these distances. However, as it was emphasized already in |2| , they are not necessary 
as precise in the limiting cases ri )2 = and ri 2 = 0. 

The motivation of our study is that in a number of dynamical problems one needs the 
bound state wave functions in the such region of variables, where one of the distances 
is much smaller then the other ones. This takes place in those processes on the bound 
electrons, which are kinematically forbidden for the free ones. For example, the high energy 
asymptotic of photoionization is expressed in terms of the two-electron function \l/(0, R, R), 
where r\ or r 2 is zero Q . The same is correct for the nonrelativistic high energy asymptotic 
of the double photoionization, and for the energy distribution of the Compton scattering at 
sufficiently small energies of the outgoing electrons. Some of the characteristics of the double 
photoionization are expressed in terms of the two-electron function with zero interelectron 
distance ty(R,R,0) 0. The straitforward way to obtain the functions F(R) and &{R) is 
to calculate them from \&(ri, r 2 , ri 2 ), that could be derived numerically. This is a rather 
complex procedure. Therefore here we shall build the approximate wave functions Fa{R) 
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and $a{R) based on the known behavior of the exact wave function only near the triple 
coalescence point R = 0. The only free parameter of our approach is the value of the wave 
functions at the coordinate origin 

N = ^(0,0,0) . (2) 

Our approach was initiated and encouraged by the large role of the proper treatment of the 

two-particle coalescence point in the earlier calculations. For example, the binding energies 

can be reproduced usually with a good accuracy by the approximate wave functions, which 

ri 

are certain combinations of exponential and polynomial factors Such presentation is 
good enough at distances of the order of the atomic size. However, it was understood long 
ago that the analytical dependence upon R is not as simple as that, and the logarithmic 
terms are presented in the expansion of the wave function near the origin j4|. Later it was 
found that, if or r i2 turn to zero, the solution of the Schrodinger equation should satisfy 
the specific Kato conditions lal. Inclusion of the logarithmic terms lol, or accounting of 

o n 

the Kato conditions |7fl or both |8|] does not influence much the energy value, but improves 
strongly the convergence of the procedure of ^(^^2,7*12) calculations. This encouraged us 
to try a rather simple approach. 

As it is shown in this paper, the approximate functions for (1) appear to be 



F A {R) = iVexp 



R 



<$> A (R) = Nexp(-2ZR) . (3) 



They have to be compared to precise or highly accurate locally correct functions Flc{R) 
and <&lc{R)- As such, we use the functions obtained by the Correlation Function Hy- 
perspherical Harmonic (CFHH) method 9]. These non- variational wave functions of the 
two-electron system bound by a light nucleus in s-state have been obtained by direct solu- 
tion of the three-body Schrodinger equation Q], without additional approximations. They 
require complicated computer codes for being employed. 

The way we construct the approximate wave functions insures that they reproduce the 
CFHH functions Flc{R) and $>lc(R) with good accuracy at sufficiently small values R. The 
question is, how long this can last, while R increases. In other words, we must calculate 
the characteristics of the processes, which are determined by F(R) and <&(R) at R being of 
the order of the size of the atom, and compare the results obtained with (3) and with the 
CFHH functions. 
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The answer is that the relative discrepancy between functions (3) and CFHH functions 
does not exceed several percents at characteristic distances ( Z -i/2) an d The same is the 
accuracy of the photoionization characteristics. 

Of course, such accuracy would not have been sufficient for the calculation of the static 
atomic characteristics, e.g. of the energy levels values. However, e.g. there was qualitative 
controversy in theoretical results on the double photoionization energy distribution until 
recent time jllj . with quantitative results differing by orders of magnitude. Thus it would 
be unjustified to run for too high accuracy in any case. On the other hand, good accuracy of 
the functions (2) prompts a basis for expansion of the CFHH functions. Since the functions 
(3) have the radial dependence of the ls-functions in the Coulomb fields with charges (Z— |) 
and 2Z, respectively, one can present the numerical CFHH functions as linear combinations 
of the functions of this field with the dominative contribution coming from Is terms. 

We build our approximate wave functions and discuss their relation to other approaches in 
Sec. 2. We analyze expansion of CFHH functions at two-particle coalescence points in series 
of the single particle eigenfunctions of Coulomb fields in Sec. 3. We consider the applications 
in Sec. 4, and summarize in Sec. 5. Atomic system of units is used through the paper. 



II. WAVE FUNCTIONS 



It is known that at small distances r± 2 C ^ 1 the solution of the Schrodinger equation 



can be presented as 



12 



3 



^(n,r 2 ,r 12 ) 



N 



1 - Z[r x + r 2 ) + -7*12 + 0(r 2 , r 2 In r) 



(4) 



with r = \Jr\ + r\. The explicit form of the quadratic terms was found in 12 1. Eq.(4) is 
consistent with more general Kato conditions [a], 

d^{r 1 ,r 2 ,r 12 ) 

= -Z^(0,r 2 ,,r 2 ); 

n=o 

= -zt>(n,o,n); 



Or 1 

d^(r 1 ,r 2 ,r 12 ) 



dr 2 

d^(ri,r 2 ,riz) 



dr 12 



r 2 =0 



r 12 =0 



-tf(n, n,0). 



which are fulfilled for the CFHH functions. Using Eq.(4) we find that at r 12 <Z 1 



F(R) = N 



Z--)R+ 



; = N (1 — 2ZR + .. 



(5) 



(6) 
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with the dots denoting the higher terms. This provides 

1 dFW = 1 
fl-o F(R) dR 2 V ; 

and 

»-o$(i2) rfi? v ; 

We require Eqs. (7) and (8) to be satisfied by our approximate functions Fa(R) and 
$ A (R) for all R. This leads to Eq.(3). 

The functions (3) correspond to a very simple physical pictures. Note that (3) look like 
the Is functions in the Coulomb fields with charges [Z — 1/2) and 2Z, respectively, that 
serve in fact as a sort of adjustable parameters. It's R dependence is one of Is electron, 
while the small probability of the three-particle coalescence is contained in the factor N 
determined by Eq.(2). We calculate the latter by using the CFHH functions. 

To characterize the quality of our approximate functions, we introduce the 

(T>\\ 1 l F A (R)-F CFHH (R) l ^A(R)-^CFHH(R) i , Q v 

y 1 (R)) = log 10 | - — 1; y 2 {R) = log 10 | — 1 (9) 

rCFHH\tt) ^CFHHy-ti) 

with the lower indices CFHH denotes the wave functions obtained in jiol ]. 

The accuracy of the functions (3) increases rapidly with the nuclear charge Z growth. 
However even for the negative ion H (Z = 1) the accuracy is rather high. At characteristic 
R ~ [Z — 1/2) -1 , and R ~ (2Z)~ l the errors of the function $^4 for H make 6%, being 
less than 1% for the function Fa- The errors increase at larger values of R. They exceed 
the value of 10% at the distances, at which the wave functions are already very small. The 
functions Vi{R) (9), describing .R-dependence of the errors are presented in Fig. 1. We 
present the results for helium [Z = 2) since most of the studies of the two-electron systems 
are carried out for this case. We give also the results for Z = 4 to illustrate Z dependence. A 
curve for H~ (Z = 1) is also presented, since this case is the most difficult for investigations. 
The dip on the graph of Fig. la is a result of the logarithmic scale, since the logarithm of 
the absolute value of the difference of the two functions goes to —00 at the points where the 
difference changes sign. The overall accuracy of the solution therefore can be inferred only 
at the values of R not too close to the dip. 

One can see that the discrepancy with CFHH functions becomes much greater at R 
becoming of the order of the size of the atom, comparing to that at smaller R. However, the 
precision is still good enough for obtaining results with the accuracy of several percents. 



The values of N, the latter being defined by Eq.(2) are presented in Table 1. At large 
Z the single-particle hydrogenlike model is expected to become increasingly true, since the 
interaction between the electrons is Z times weaker than their interaction with the nucleus. 
Hence, in the limit Z ^> 1 

N = N c = — . (10) 

7T 

The results presented in Table 1 illustrate this tendency. As expected, deviations from the 
limiting law (10) are of the order Z~ l . The actual results are smaller than predicted by (10) 
since the latter does not include the electron repulsion, which diminishes this value. 
Of course, there are numerous simple approximate wave functions of the type 

^(''l, r 2, 7"i 2 ) = c(exp(-ar 1 - br 2 ) + exp(-ar 2 - br x )), 

which are build in order to calculate the ground state energy values fjj , thus approximating 
the solutions of the Schrodinger equation at ri 2 of the order Z~ x (in the case of H~ they 
had to reproduce also the very existence of the bound state). Technically, they turn to the 
single-exponential forms at r\ = r 2 = R and do not depend on 7*12. These functions can 
be compared to our functions ®(R) defined by Eq.(3). But they do not approximate the 
locally correct CFHH functions &a(R), and, following |2j, are not supposed to. We illustrate 
this statement by presenting in Fig.2 the CFHH function $>{R), our function (3) and the 
screened Coulomb wave function <& S (R) = —exp(—2aR) with a = 27/16 for helium 

In [14I the function F(R) for H~ , He and Li + was approximated by a hydrogenlike 
function with the effective charge Z e tt treated as a variational parameter. The values of 

r 

Z e ff for Z — 1,2,3 have been found to be 0.58, 1.53 and 2.52, correspondingly. In |15j the 
function F(R) for the ion H~ have been analyzed at large distances. We do not claim our 
functions to be true in this i?-region, which is not essential for us since of prime importance 
is the R domain within the atomic radius. 



III. EXPANSION OF CFHH FUNCTIONS IN SERIES OF THE COULOMB 
FIELD EIGENFUNCTIONS 

The R dependence of the approximate wave functions Fa(R) and <&a(R) (3) is the same 
as that of Is functions in the Coulomb fields of the nuclei with the charges Z\ = Z — 1/2 and 
Z 2 = 2Z, respectively. The high precision of these functions suggests that the eigenfunctions 
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of the Schrodinger equations in these fields compose convenient series for expansion of the 
CFHH functions F(R) and $(R). 

Introducing the common notation X(R) for the functions F(R) and we present the 

oo oo 

normalized functions X N (R) = -^X(R) with C x = J R 2 X 2 (R)dR. Thus J R 2 X 2 ! (R)dR = 

C x o 

1. 

In the expansions over the complete sets of some eigenfunctions, X N (R) can be presented 

as: 

F N (R) = £a^(i?); <S> N (R) = £& <¥ > i (i2) , (11) 
with denoting the sum over the states of discrete spectrum and integration over continuum 



oo oo 



a,i = J R 2 F N (R)f*(R)dR; bt = J R 2 $ N (R)<p*(R)dR . (12) 



For fi(R) and <fi(R) normalized to one, it is 

£«? = £** = !■ (13) 

Choosing the solutions of the Schrodinger equations in the Coulomb fields with the charges 
Zi = Z — 1/2 and Z 2 = 2Z as the functions fi(R) and (fi(R) respectively, we find the values 
au and bi s presented in Table 2. For atomic helium ai s = 0.9997, b\ s = 0.998. High accuracy 
of the functions (3) corresponds to domination of the terms a\ s and b\ s in the sums (13). 

The precision of calculations can be improved by adding the contributions of the higher 
states according to Eq.(12). Of course, in our case only the s-states are involved. For 
example, a2 S = —0.02, 6 2 s — —0.05 in the case of atomic helium. The results for the other 
values of Z are presented in Table 2. This procedure enables to achieve any desired accuracy, 
controlled by Eq. (13). 

IV. EXAMPLES OF APPLICATION 

As we said above, one of the possible application of the functions (3) is the high energy 
photoionization processes. Let us start with the single photoionization. The high energy 
nonrelativistic asymptotic for the K-shell ionization cross section can be written as |l( 

2 n / 2 7ie 2 Z 2 C 2 
3mcuj 7 / 2 
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where m is electron mass ans c is the speed of light. The properties of the ionized states 
contained in the factor 



oo 



C = J R 2 F{R)ip K {R)dR . (15) 

Q 

Here F(R) is determined by Eq.(l), while ipxiR) is the single-particle function of the K- 
electron in residual ion. In our case %Pk{R) is just the Is function of the Coulomb field with 
the charge Z. 

In the single-particle approximation C is simply the value of the single-particle wave 
function at the coordinate origin. To illustrate the quality of the functions (3) we compare 
the results for the factor C calculated by using the CFHH functions and the functions (3). 
In the latter case we find an analytic expression 

2NZ 3 / 2 

°~ v^F(2Z- 1/2)3 > ( 16 ) 
providing C = 0.102 for the case of atomic helium. The numerical calculations with the 
CFHH functions give C = 0.103 in this case. Hence, employing the approximate function 
(3) leads to the error of 1%. Earlier the authors of |l6( found that the value of C ob- 
tained by using the Hylleraas-type variational function is well approximated by employing 
a hydrogenlike function with Z e ff = Z — 0.53. 

Now let us turn to the case of the double photoionization. The shape of the spectrum 
curve of the double photoionization changes with the photon energy growth. The mecha- 
nisms which cause these changes are explained in While the photon energy uj is smaller 
than certain value ui, the energy distribution approaches its minimum at the central point, 
with the equal energies of the outgoing electrons, e^2, i-e. S\ = e<i- There is a peak at the 
central point at uj > u>%, which splits into two at uj > uj 2 . Thus, there is a local minimum at 

E\ = E2 at LO > U02- 

The values of were obtained in 17( by using the CFHH functions. We shall not 
repeat derivation of the corresponding equations here. We rather explain their origin and 
put them down, in order to illustrate, how the functions (3) enable to obtain approximate 
solutions. 

The values of ui^ can be presented as solutions of the following equation, which involves 
the functions F(R) and $(R) Q: 

X/i = uj 9/2 A(uj) , (17) 



with A being a certain numerical coefficient, and 

oo 

H = J dr\F(r)\ 2 , (18) 
o 

while the function A depends on w in a more complicated way: 

+1 

A(u) = J dtt 2 (l - 2t 2 )D(uj 2 t 2 ) (19) 



with 



D(q 2 ) = \ I S ^l^( r )r 2 dr\ 2 . (20) 



Employing the exact CFHH functions requires tedious computations. However, using 
the approximate wave functions (3) one can obtain analytical expressions for both left-hand 
side and right-hand side of Eq.(17). Putting F(r) = Fa(t) and $(r) = $a(^) we obtain 
H = 2^rY > while 

1 /6a 6 + 13a 4 + 2a 2 + 3 1 - 2a 2 1\ 

Aiuo) = — —. - — 1 — - — arctan- (21) 

v 1 uj g V 6a 2 (a 2 + l) 3 2a 3 a] v ' 

with a = IZjuj. 

The values of uip obtained by using the CFHH functions and the functions (3) are 
presented in Table 3. One can see that the discrepancy between two sets of results drops 
rapidly with Z growth. Being 22% for H~ and 9% for helium, it becomes 4% for Z — 4. 



V. SUMMARY 

We build very simple analytical presentations (3), for the wave functions F(R) and $(i?) 
describing ground states of two-electron systems bound by the Coulomb field of the nucleus 
in the space regions ri j2 = and r±2 — 0. The presentation is based on the behavior of the 
exact solution of the Schrodinger equation near the three-particle coalescence singularity. 
Comparing our functions (3) to the locally correct CFHH functions for the ion H~ , atomic 
helium and light heliumlike ions (relativistic corrections, which are of the order (Z/137) 2 
are not included), we found good agreement in a large interval of the values of the R. As 
is evident the precision of the approximate functions increases with the nuclear charge Z 
growth. 
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We show that the solutions of the single-particle Schrddinger equations in the Coulomb 
fields with the charges Z\ = Z — 1/2 and Z 2 = 2Z provide natural basis for expansion of the 
functions F(R) and &(R) with the domination of Is terms. The latter tendency increases 
with Z. The approach is more precise for the function F(R), then for <&(R). 

Examples, presented in Sec. 4 show that even for the lightest heliumlike systems such as 
H and He the wave functions (3) can be used at least for the estimation of the physical 
parameters. 

The high precision of such a simple approximation that properly treats singularities in the 
wave function is in agreement with the conventional believe that the singularities determine 
such important atomic characteristics as high-energy photoionization cross sections. 
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TABLE I: The value F(0) = $(0) = N for several values of Z. The ratio r = N/N c with N c 
defined by Eq. (10) illustrates the convergence to the high Z limit. 



z 


1 


2 


3 


4 


5 


6 


N 
r-N 


0.071 
0.22 


1.37 
0.61 


5.77 
0.67 


15.2 
0.74 


31.6 
0.79 


56.8 
0.83 



TABLE II: The coefficients of the two lowest terms of expansions (11) of the CFHH functions in 
terms of the Coulomb functions. The coefficients of the next terms are limited by the conditions 
|a.;| < a, \b,i\ < b, while the values of a = (1 — a\ s — a\ s ) l l 2 and b = (1 — b\ s — b\^) x l 2 are presented 
in the two bottom lines. 



z 


1 


2 


3 


4 


au 


0.98482 


0.99970 


0.99991 


0.99996 


hs 


0.99067 


0.99807 


0.99918 


0.99955 


a.2s 


-0.144 


-0.020 


-0.010 


-0.007 




-0.108 


-0.046 


-0.030 


-0.022 


a 


0.097 


0.015 


0.008 


0.005 


b 


0.082 


0.041 


0.028 


0.021 



TABLE III: The values of oj\ and L02 (Sect. 4) in keV for the ground states of the lightest heliumlike 
systems, calculated by using the CFHH functions [lfj and the functions (3). 



z 


1 


2 


3 


4 


uji - this work 


0.67 


2.11 


3.92 


6.14 


wi - [10] 


0.55 


1.93 


3.70 


5.89 


UJ2 - this work 


4.86 


9.71 


14.5 


19.3 


UJ2 ~ \M 


3.97 


8.89 


13.7 


18.5 
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FIG. 2: The exact and the approximate helium wave functions at the electron-electron coalescence 
line. 
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